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Abstract 

The kinematic reconstruction of neutral current high events at HERA is 
discussed in detail using as an example the recently published events of the HI 
and ZEUS collaborations at > 15000 GeV^ and M > 180 GeV, which are 
more numerous than expected from Standard Model predictions. Taking into 
account the complete information of these events, the mass reconstruction is 
improved and the difference between the average mass of the samples of the two 
experiments is reduced from 26ibl0 GeV to 17ib7 GeV, but remains different 
enough to render unlikely an interpretation of the excess observed by the two 
collaborations as originating from the decay of a single narrow resonance. 



1 Introduction 



At the electron-proton {ep) HERA collider, the study of DIS is performed in an unique 
and optimal way up to values of the squared momentum transfer of 10^ GeV^. 
Recently, the two HERA collaborations HI and ZEUS have reported Q an excess 
of events at very high (defined in the following as > 15000 GeV^) compared to 
standard model expectations. This observation has triggered an important activity on 
"Beyond the Standard Model" theories which might explain the effect 0. The most 
favoured solution involves the production of a new resonance, which after its decay is 
kinematically similar to a standard deep inelastic scattering (DIS) event. 

Indeed, in the naive quark-parton model (QPM), a parton carrying a fraction x 
of the momentum of the proton scatters elastically against an electron and in the 
Standard Model, this interaction is representing a t-channel scattering occuring via the 
exchange of a gauge boson such as a photon or a Z°. Beyond the Standard Model, the 
naive QPM picture can represent the formation of an s-channel resonance (generically 
called "leptoquark" ) subsequently followed by a two-body decay. While in the first 
interpretation the Bjorken x variable is one of the two relevant variables to characterize 
the scattering, in the second case the invariant mass M of the system formed is the 
physical quantity of interest. This mass is related to x in the naive QPM by M = 
^/xs, s being the squared center of mass energy. In "real" interactions, the quantum 
chromodynamics (QCD) effects do not modify much this picture at large Q^, since it 
is characterized in DIS by the dominance of one electron -|- one partonic jet -|- one 
remnant jet in the final state, and is thus similar to the one in which a leptoquark is 
produced and decays in 2 bodies. So, when looking at high Q"^ deep inelastic scattering 
(DIS) in an inclusive way it is possible to distinguish between these two scenarios by 
comparing the event rates in different x or M intervals: the smooth evolution of the 
DIS cross-section as a function of x can be opposed to the appearance of a sharp peak 
in the invariant mass distribution, characteristic of a leptoquark. 

The HERA effect is complicated by the fact that the invariant mass distributions 
of the event samples of the two collaborations differ: at a mass greater than 180 
GeV, the 7 HI events appear clustered around a mass of 200 GeV, while the 5 ZEUS 
events are broadly distributed between 191 and 253 GeV. It is the subject of this 
paper to understand how significant is this difference and to provide a combined mass 
spectrum using the complete available information of these events in order to give an 
interpretation of the effect. 

At HERA the kinematic reconstruction does not need to rely on the scattered lepton 
only, since the most important part of the hadronic system is visible in the almost An 
detectors HI and ZEUS. This redundancy allows for an experimental control of the 
systematic errors and the radiative corrections, hence to determine in a more precise 
way the usual DIS kinematic variables x, y and which are defined as: 

x = gV(2P.g) y = {P.q)/{P.k) = -{k - k'f = -q^ = xys 

P, k being the 4-vectors of the incident proton and lepton, k' the scattered lepton one. 
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In this report we briefly review in section 2 the methods used at HERA to determine 
the kinematics of the high events. In section 3 we characterize the high region 
and introduce a new method which optimizes the kinematic reconstruction by deter- 
mining the measurement errors on an event by event basis. Section 4 is devoted to the 
detailed study of the very high HERA events and discusses the mass distributions 
obtained. 



2 Kinematic Reconstruction at HERA 

In order to introduce the kinematic methods used at HERA in the high region 
let us start by a few definitions. The initial electron and proton (beam) energies 
are labeled Eo and Pq- The energy and polar angleQ of the scattered electron (or 
positron) are E and 9. After identification of the scattered electron, we can reconstruct 
the following independent hadronic quantities: S, obtained as the sum of the scalar 
quantities Eh — Pz,h associated to each particle belonging to the hadronic final state, 
PT,h as its total transverse momentum and define the hadronic inclusive angle 7: 

^ = Y,iEh-pz,h) PT,h = J{Y.P^,hy + iY.Py,hy tan^ = ^ (1) 
h \l h h ^ PT,h 

Eh,Px,h,Py,h,Pz,h are the four-momentum vector components of every energy deposit in 
the calorimeter originating from any hadronic final state particle. The corresponding 
quantities for the scattered electron are 

6 Yj 

Eg = -E (1 — cos ^) pTe = Esm6 i.e. tan - = — — (2) 

2 pT,e 

Out of these variables, it is possible to write the four methods which have been used 
to determine the invariant mass of the very high events (the formulae are given in 
the appendix). 

- The electron only method (e) based on E and 6. 

- The hadrons only method (h) based on E and pT,h @]- 

- The double angle method (DA) based on 6 and 7 0. 

- The Sigma method (S) based on E, 6 and S 

The e method is precise at high y but becomes less precise in x,y 01 M at low y. 
The h method is the only one available for charged current events, but is less precise 
than the three others in neutral currents, so it will be ignored in the following. The 
DA method is precise at high both at high and low y. The S method is similar to 
the e method at high y and to the DA at low y apart from which is less precise. 
However it is insensitive to initial state QED radiation (ISR) in M, y and Q^. See 
for a detailed discussion of the properties of these methods and ||^ for their behaviour 
in presence of ISR. 

^The positive z axis is defined at HERA as the incident proton beam direction. 
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It should be stressed that the quahty of any kinematic reconstruction depends 
mainly on the precision achieved on the observables, such as E, 9, S or 7. It is one 
of the major tasks of the experiments to obtain precise and unbiased quantities, and 
several techniques have been developped in the HERA structure function analyses to 
achieve these goals. We refer to the original publications for a description of these 
studies 0, which discuss for instance how the electron energy is calibrated, or how the 
hadronic final state is measured at low y, when the hadronic jets are in the vicinity of 
the beam pipe. 

3 Treatment at High 

At high (defined in the following as > 2500 GeV^) the kinematic reconstruction 
is in general more precise than at low and indeed the differences between the results 
obtained with the methods seen above are small. However since the number of events 
drops rapidly with increasing it is still crucial to optimize the reconstruction by 
making a full use of all the observables of these events. 

The improvements on the reconstruction at high come mainly from the better 
measurement of the hadronic final state. This is due to the fact that i) the individual 
hadron energy is on average greater; ii) the losses in the beam pipe and in the material 
in front of the calorimeter are in proportion smaller than at low Q^, i.e. S, pT,h and 
7 are less affected by these losses; iii) the hadronic final state displays more often a 
(single) collimated jet configuration. This is due both to kinematics (events are in 
average at higher x), and to the smaller gluon radiation and power corrections {as is 
smaller) . 

All these characteristics enable a precise measurement of the hadronic angle, and 
thus are particularly favourable to the DA method which indeed becomes more accurate 
with increasing Q^. This method has however the drawback to be very sensitive to 
ISR and it is thus difficult to apply to a small number of events. For instance, the 
ISR of a 2.75 GeV photon, if not taken into account produces a shift of 10% on the 
reconstructed mass and 20% on Q^. In order to overcome this drawback we introduce 
a new method which, by making use of the precision on 7, estimates from the data 
on an event by event basis the shift in the measurement of E and S, denoted in the 

following SX/X = {Xtrue- Xreconstructed) / Xreconstructed) with X = E OI ThcSC shiftS 

allow i) the E method to be corrected; ii) the presence of an undetected initial state 
photon (down to an energy of about 2 GeV) to be identified with high efficiency, and 
therefore to correct for it. 

3.1 The uj Method 

Let us assume that there is no ISR (the specific treatment of ISR is discussed in the 
next section, although the two steps are not separated in the procedure), and that 6 
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and 7 are precisely measured (implying = 6pT,h/pT,h)- From energy momentum 

conservation the two following equations can be derived: 



[i-ye) ^ + yh ^ = ye- Vh (3) 

6E 

-pT,e + PT,h = PT,e " PT,h (4) 

Under these assumptions 6E/E and are determined on an event by event basis. 

Varying 6 and 7 within their typical errors^ (5 and 40 mrad (r.m.s.) respectively, at 
high Q"^) we can obtain the errors arising from the angular measurements on 6E/E and 
(5E/E. A more direct way to see the uncertainties arising from this determination is 
illustrated in fig.l, in which the error on E and S reconstructed using equations 3 and 
4 are compared to their "true" value, which is known for the reconstructed simulated 
events. For a given event, the error on {6E/E)rec, which is defined as A{6E/E) = 
{SE/E)rec — {SE /E)true, is of the same order as the r.m.s. (denoted < 6E/E >) of the 
SE/E distribution obtained from a large sample (fig. lb), showing that its use will not 
bring an improvement on an event by event basis. On the other side, the relative error 
on ((5S/S)rec is smaller, about 30% of < > , and diminishes at high y (compare 

fig.le to If, and lb to Ic) implying that correcting E will provide a better kinematic 
measurement, in particular at high y. 




(6I/E)-rec 



&(6T/T)/<6T/L> 



A(6I/I)/<(5I/I> 



Figure 1: a) shows the correlation between true reconstruction errors on E and those 
obtained from eqs. 3 and 4, on a full simulation^ in the HI detector of high events; 
b) shows the projection of this correlation; c) shows the same projection restricted to 



y > 0.4. We use A(^) = (^)rec - (^)ir«e, and < ^ >= r.m.s. of the ^ distribution 

S 



3%. d,e,f) show the same plots but for E instead of E, with < ^ > = 10%. 



^At very high Q^, the 7 resolution improves to 30 mrad. We use here the r.m.s., rather than the 
standard deviation obtained from a gaussian fit to the distribution, to take into account non-gaussian 
tails which contribute in the propagation of the systematic errors. 
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The oj kinematic variables are thus derived from the S ones by including the effect 
of 5S/S, i.e. 

" S + 5S + Se " 1 - ^ 

The comparison of the uj and S reconstruction of M, y and for high events fully 
simulated^ in the HI detector is shown in fig.2. The improvement obtained by the 
recalibration of E is clearly visible, allowing this method to be compared favourably 
to the original ones. This comparison is shown in fig. 3 for high events (in all the 




Figure 2: The improvements on the kinematic reconstruction at high , when using 
the uj method compared to the S method, are shown for the mass M (2a), y{2b) and 
Q\2c). 

high section, an E — cut against hard initial state radiation is applied, see next 
section) both at high y and at low y. At high y {y^ > 0.4), the e and the u methods are 
comparable in M and y and slightly better than the DA one (fig.3a,b). In Q^, the DA 
displays the narrowest peak both at high and low y (fig.3c,f), but its high sensitivity 
to ISR induces tails in the distribution which renders its r.m.s. larger than the u and 
e ones. At low y [y^ < 0.25) the e method has a relatively poor resolution in M, much 
worse than the uo one which is slightly better than the DA one. In conclusion, for the 
high events the uj method is similar or slightly better than the e and DA method. 



3.2 Treatment of QED Initial State Radiation 

The influence of ISR is taken into account in the simulation programs used by the 
HERA experiments, but on a small number of events it is difficult to control the 
migration due to an unseen ISR photon when using the e and DA methods. Both HI 

^In figures 1 and 2, the radiative processes of the DJANGO program which is used to generate 
the events have been turned off. The complete simulation, based on DJANGO, ARIADNE ||ll| and 
GEANT |]l2| as used in the HI simulation program, is further described in [Q. 



5 



2800 



21 00 



1 400 



700 







(J 

e 

DA 



1 500 



1 125 — 



750 — 



375 — 



1 5200 r 

mean rms _ (b) 

1,30 0.04 

0,99 0-04 2400 — 

1.00 0.05 Z 

1600 — 



800 — 



neon rms 

1.01 0.05 

1.02 0-05 
1,01 0.05 



2400 



1 800 



(c) 



mean rms 

1.00 0,07 
0,99 0,06 

1.01 0.08 



1200 — 



600 — 



0.6 0.8 1 1.2 1.4 



0.8 1 '.2 1.4 

Yrec/' Ygen 








0.6 0.8 1 1.2 1.4 



mean rms 




1 000 



750 — 



500 — 




2800 



(e) 



250 — 



2100 — 



1400 — 



700 — 



mean rms 

1.00 0.06 

0.99 0.06 

1,31 0,07 



J; 



0.6 0.8 1 1.2 1.4 0.6 0.8 1 ',2 1,4 



0.6 0.; 



1 1.2 1.4 

Q rec/ Q ger 



Figure 3: Comparison on high events, at high y (ye > OA) of the uo, DA and e 
reconstruction of the mass (3a), y (3b) and (3c). The same comparisons are made 
at low y (ye < 0.25) in fig. 3d,e,f. 



and ZEUS use an experimental cut on the total E—p^ of the event (T^he) or equivalently 
on the normalized E — called hereafter ahef\ and defined as 

_ S/ie _ S + Se _ , . 

f^/^e = ^ = =(Th + (Te (6) 

The E—pz cuts used in HI and ZEUS correspond to good approximation to ahe > 0.75 
and prevent an ISR photon from carrying away more than about 25% of the incident 
electron energy. Such a photon can induce a large shift to the reconstructed mass M of 



1 - ^(0.75?/e)/(l/e - 0.25) for Me (i.e. -18% at ye = 0.75) and of 25% on Mda at any 
y. On Ms the shift is zero since the method is independent of colinear ISR for M, y 
and Q^. A more stringent cut on ate is difficult to implement due to the experimental 
resolutions on yt and ye- 

The u method allows ISR to be identified with high efficiency for a^e as high as 
~ 0.93. It makes use of the simple fact that depending on the origin of the ate shift 
(which comes mainly either from ISR or from hadronic miscalibration) , the errors 

^The normalized E-p^ for the electron and for the hadronic final state satisfy: Ue = 1 — ye, <Jh — Uh- 
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obtained with the u method have a completely different pattern: for example if the 
observed ahe is 1— 2, equations 3 and 4 will give for the ISR case 6E/E = z and 
(5S/Z1 = z, assuming negligible detector smearing, while for the hadronic miscalibration 
case, assuming negligible error on the electron energy, we will get 6E/E = and 

= z/uh- These two different types of pattern allow ISR to be identified even 
in the presence of detector smearing, as studied on a complete simulation in the HI 
detector, with typically 85% efficiency at high for photon energies greater than 2 
GeV. The exact conditions for ISR identification of a 2 GeV photons is 6E/E and 

> 7.3%. However, to take into account the detector smearing we use instead: 

— > 5% and — > 5% and — + — > 15% (7) 
E S E Yj 

These can be slightly varied depending on the ISR identification efficiency requested. 
These conditions will "misidentify" high non-radiative events in less than 3% of the 
cases in the HI detector, implying that in the total sample of ISR identified events, the 
fraction of genuine ISR events is higher than the fraction of non-radiative events. This 
fraction increases with increasing since the hadronic angle becomes more precise, 
i.e. the ISR identification becomes more efficient. In case of ISR identification, the 
equations 3 and 4 are solved again after recalculating i/e and i/h using ahe ■ 2-Eo instead 
of 2Eo. The improvement obtained can be visualized in fig. 4 which shows the same 
distributions as in fig. 3 but on the sample of DIS high events having radiated an 
ISR photon of 2 GeV or more. The dramatic improvement underlines the importance 
of controlling experimentally the radiation effect. Furthermore, the soft ISR photons 
which cannot be identified below 2 GeV have a small effect on (< 2%) but still a 
sizeable one on Mda, up to about 7%, or 4% on Me. 




Figure 4: Comparison of the uj,e and DA reconstruction of the mass (4a), y (4b) and 
Q"^ (4c) for high Q'^ radiative events in which El^^^ > 2 GeV (and ate > 0.75). Note 
the similar behaviour ofy^ and yoA, since yoA is ISR independent. In M and the 
improvements brought by the u reconstruction method are remarkable. 
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4 The Very High Events at HERA 



Let us remind some characteristics of the excess of events at high observed by the HI 
and ZEUS collaborations in the data collected from 1994 to 1996 |l], g. At > 15000 
GeV^, 24 events have been observed for an expectation of 13.4±1.0, representing a 



probability of 0.0074 [13]. Furthermore some of these events exhibit the additional 
characteristics either to cluster around a rather precise value of the invariant mass of 
the produced system (Me = 200 ± 12.5 GeV, for 7 events out of the 12 of HI, while 
1.0 ± 0.2 are expected), or to lie at very high invariant mass (M/j^i > 223 GeV and 
Uda > 0.25, for 4 events out of the 12 of ZEUS, while 0.9 ±0.2 are expected). Actually 



the messages of these two observations could be different as already pointed out |1^, |T5 
and we will try here to explore further this point with the use of our new kinematic 
tool. 

Among the obvious possible reasons for the difference are the influence of ISR 
which can strongly distort a distribution based on a small number of events, or the 
effect of a speciflc miscalibration. The u method having been devised to respond to 
these problems, we can now study in more details the 11 extreme events mentioned 
above for which the complete kinematic information has been published by the two 
collaborations. We add to this sample 2 events^ from ZEUS in order to include all 
events which satisfy the condition > 15000 GeV^ and M> 180 GeV, since we want 
to test the signiflcance of the HI event clustering at masses around 200 GeV over the 
complete very high HERA sample. 



4.1 Kinematic Properties 

In table 1 are given the M, y and of the 13 events at very high and high mass 
reconstructed with the four different methods discussed above. The values of M, y and 
are identical to those published for the e and DA methods. The S values have 
been released recently by HI |T6[ and can be computed straightforwardly for the ZEUS 
events. The errors are also reproduced exactly except for those of the e method to 
which we conservatively added quadratically an error of 1.5% (i.e. half of the absolute 
energy scale uncertainty of 3% given by the two experiments) to account for potential 
miscalibrations between the different regions of the detectors. As we will see below the 
overall energy calibration can be checked to be well under control. 

For the errors not published by the collaborations {u, S for HI and ZEUS, and DA 
for HI) we used the following prescription, which was checked to be consistent with 
the other published errors: They are obtained here using a possible error of ±5 mrad 
and ±30 mrad for the electron and hadronic angle respectively. For the error on the 

^For one of these 2 events (called Z-6 in table 1) only the DA variables are available (from fig.l 
of rcf. Q). The additional kinematic properties were deduced from the DA values using the average 
E — pz and systematic shifts between e and DA variables of the 5 other ZEUS published events. This 
assumption has a negligible influence on the conclusions drawn below. 
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electron (hadronic) energy we used ±3% (±4%) for HI and ±5% (±4%) for ZEUS, 
which includes for the electron case both resolution effects, dead material corrections 
and the 1.5% overall error just mentioned. These two values for the electron energy 
error have been checked on the data (sec below). Additional sizeable errors as published 
by the collaborations, due for instance to special energy corrections, have been taken 
into account in two events (H~3,Z-1). 

For the HI events the errors of an event are similar for each method, with a slight 
advantage for the e method at y above 0.5. For the ZEUS events the DA errors are the 
smallest, since the error on the electron energy which is somewhat larger than in the 
HI case has an influence on the 3 other methods. This also explains why the ZEUS 
collaboration chooses the DA method rather than the e method favoured by HI. 

Also mentioned is the result of the uo ISR idcntiflcation, and we can see that 2 
events (H-5 and Z-4) are classified as radiative, and have been corrected accordingly. 
These 2 events indeed show the characteristics of a radiative event, i.e. there exists a 
value ol z = E^^^/Eo for which the following equations, which hold exactly if there 
were no detector effects, are verified in a good approximation: 

5E/E ~ z ~ z (8) 

Ql/Ql ^1-z QIIQIa ^ 1 - ^ (9) 

Ms/M^A '^l-Z Xj,/XDA ^ (1 - zf (10) 

Vda ^Vt. ahe^l- z (11) 

Indeed in the case of H-5 a photon with an energy (2.9 GcV) consistent with the 
z derived from the previous equations is observed in the special photon calorimeter 
located close to the beam pipe and designed to measure photons emitted coUinearly 
to the incident electron, such as ISR photons. The ISR effect on the mass is sizeable, 
since M^a/M^ = 1.08 for H-5 and 1.11 for Z-4. Both values are consistent with the 
determined 6E/E and of these events (sec tab.l). 

The reconstructed kinematics are consistent between the different methods except 
for the 2 radiative events and for the events displaying a significant difference between 
Me and Mda (events H-7,Z-l,Z-5), which indeed have the largest 5E/E and/or 5E/E. 

On these small samples of events, and after ISR corrections, the averages of the abso- 
lute value of the error on the hadronic energy are similar for HI and ZEUS (< |^| >= 
6% compared to 5%), while on the electron energy the HI average error is smaller: 
< |M| >^3% compared to 5% for ZEUS. 

The uncertainty on the electromagnetic absolute energy scale can be estimated from 
the data using < ^ > (the 2 radiative events are excluded from these means) which 
gives +1.1% ± 1.5% for HI and +2.1% ± 3% for ZEUS, both values in good agreement 
with the value of 3% given by the collaborations. The uncertainty on the hadronic 
absolute energy is obtained similarly and gives —5.4% ± 3% for HI and +3.3% ± 4% 
for ZEUS, also in acceptable agreement with the values quoted by the collaborations. 

After all these consistency checks, the use of the u method in the two samples will 
now allow a consistent mass distribution to be derived from the very high events. 
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.299 
.042 


15.9 
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16.6 
0.5 


16.2 
0.9 


15.8 
1.6 


Z 4 

0.92* 


+.12* 
+.07* 
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9 


253 
6 
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12 
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10 


.728 
.014 


.721 
.008 


.752 
.021 


.731 
.021 


37.5 
2.5 


46.1 
1.6 


41.0 
3.1 


37.8 
3.4 


Z 5 

0.94 


+.10 
-.03 


207 
13 


232 
10 


200 
15 


206 
15 


.305 
.024 


.285 
.017 


.350 
.033 


.310 
.044 


13.1 
1.2 


15.4 
0.4 


14.0 
0.7 


13.2 
1.4 


Z 6 

0.95 


+.03 
+.07 


185 
12 


191 
11 


186 
12 


183 
12 


.608 
.023 


.592 
.028 


.610 
.054 


.591 
.054 


20.7 
2.1 


21.6 
1.6 


21.0 
1.5 


19.8 
3.2 



Table 1: Kinematic properties of the 13 events observed by the HI and ZEUS collab- 
orations which have, at least in one method, > 15000 GeV^ and M > 180 GeV. 
For event Z-6, only the DA values are accurate, the values of the other methods being 
extrapolated. The values of M, y and are identical to those of the original papers 
(e, DA, E for HI ^, |T§, e, DA for ZEUS The new values of the table are for the S 
(ZEUS) and for the uo (Hl+ZEUS), and for the non available errors (see text for more 
details). Also given are a^ej^cind which allow to quantify the size of the errors and 
to tag the presence of initial state radiation and to correct for it (marked with a 
After ISR corrections, the (ahe,^,^) values are for event H-5: (1.00, .00, .00) and for 
event Z-4: (1.00,+.04,-.01). The masses are given in GeV, the Q'^ in 10^ GeV^. 
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4.2 Mass Distribution 



In the following we will not consider any more the event Z-5, since it survives the 
cut only for the DA method. It might be a radiative event which could not be identified 
due to a large smearing in the hadronic energy. In any case (radiative or not) it would 
not survive the cut, which is one of the 2 conditions the final sample must satisfy. 
Note "en passant" that its mass is reduced from M£)yi=232 GeV to M^=207 GeV. The 
mass distributions of the remaining 12 events, obtained with the e, DA, S and uj, are 
displayed in fig.5a,b,d,e. 
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Figure 5: Mass distributions obtained with the e (5a), DA (5b), S (5d) and u (5e) 
reconstruction methods. The main histogram represents the mass distribution of the 
HI and ZEUS events at > 15000 GeV^ and M > 180 GeV^. The dashed histogram 
represents the contribution of the ZEUS events only. The grey histogram represents 
the prediction from Standard Model DIS, obtained from a simulation normalized to the 
summed luminosity of HI and ZEUS. Fig. 5c displays the u mass reconstruction of a 
200 GeV narrow resonance generated with the LEGO program / [T^/ and fully simulated 
in the HI detector. Fig.5f is obtained from fig.5e by filling one gaussian distribution 
per event, with mean and width given by its co-mass and its error (cf table 1). 
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The event Z-4 has the highest mass at 253 ± 6 GeV. It is identified as a radiative 
event and its mass is reduced by the uj determination to 226 ± 5 GeV, thereby reducing 
the scattering appearance of the ZEUS events. Note that the error quoted on the DA 
mass is much smaller than the effect due to radiation. 

For the 5 ZEUS events the weighted average mass is decreased from M^^=226±9 
GeV to M''J'3=21Q±7 GeV. Actually with the u determination 4 out of the 5 ZEUS 
events lie between 220 and 228 GeV, the S*'^ one (Z-6) being at 185 GeV, i.e. at a 
lower mass than any of the 7 HI events (see fig.5b). 

The average mass0 of the 7 HI events is essentially not dependent on the method 
used: M^''^'=199±2.5 GeV, M^''5=200±2.6 GeV, M^''l=201±5 GeV. The 7 HI events 
remain clustered between 188 and 211 GeV. Thus the HI and ZEUS samples are con- 
centrated at significantly different mass values and this splitting cannot be accounted 
for either by ISR or by detector effects. 

The small number of events involved prevents a definite interpretation of this effect. 
However the fact that no event among the 5 of ZEUS is found in the bin where the 7 HI 
events are measured suggests that this specific accumulation is a statistical fiuctuation. 
In fig. 5c is shown the reconstruction of a narrow resonance generated at M=200 GeV 
with the LEGO [|T7| program in the HI detector using the uj methodQ. The width of 



the distribution which include experimental and QED/QCD radiation effects, cannot 
accomodate the tails of the measured experimental distributions (fig.5a,b,d,e). 

If we make an ideogram from the histogram 5e i.e. if we apply gaussian smearing to 
each of the 12 events using their uj mass as a mean, and their error as r.m.s., we obtain 
fig.5f, which also shows the incompatibility with fig. 5c. To reconcile the 2 "peaks" 
visible in fig.5f would require to miscalibrate uniformly the electron energy of the HI 
events at least by +6% and the ZEUS events by —6%, values completely incompatible 
with the absolute scale uncertainties found in the previous section. Note however, that 
the HI events alone support the narrow resonance hypothesis. 

Since after the uj kinematic treatment none of the 12 events migrates outside the 
very high and mass region, we confirm that the visible excess published by the 
collaborations at very high is not due to detector (calibration) or radiation (ISR) 
effects. 



^ The fact that the average masses are shghtly different from the original pubhcations is due to 
the difference in the errors which weights the events in a different way. This difference is irrelevant in 
the current discussion. 

^AU the methods give similar distributions on these events, except for the DA which has a larger 
r.m.s. due to radiative tails. In units of GeV the (mean; r.m.s.) are for the e, DA and E respectively: 
(195.5; 7.5), (195.3; 7.4), (195.2; 9.6), (195.3; 7.8). The bias of about 5 GeV can be removed by taking 
into account the mass of the jet, but is present when using inclusive methods. 
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5 Conclusion 



We have introduced a new reconstruction method ("a;") which allows the kinematic 
variables of the high Cf" events to be determined in a more precise way. It uses the 
kinematic constraints on an event by event basis to calibrate the hadronic energy and 
to identify and correct for the presence of QED initial state radiation. 

This method has been applied to the 12 HERA events observed by the HI and 
ZEUS collaborations (for an expectation of about 5) at > 15000 GeV^ and M> 180 
GeV. The accumulation of the 7 HI events around 200±12.5 GeV is confirmed. The 
average mass of the 5 ZEUS events is decreased from the double-angle value of 226±9 
GeV to 216±7 GeV. However none of the ZEUS events enter the HI accumulation 
region. Taking the estimation of the systematic errors at face value, this suggests that, 
if the observed excess at very high is due to physics beyond the Standard Model, 
it is unlikely to be explained by the decay of a single narrow resonance such as a 
leptoquark. The increase of the integrated luminosity in the current and in the future 
years at HERA will clarify the origin of this effect. 
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Appendix: y and formulae (M = ^jQ^/y) 
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